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Abstract 

Results are reported for an investigation into the methods 
by which energy transfer is calculated in the Direct Simula- 
tion Monte Carlo method. Description is made of a recently 
developed energy exchange model that deals with the trans- 
lational and rotational modes. A new model for simulating 
the transfer of energy between the translational and vibra- 
tional modes is also explained. This model allows the vi- 
brational relaxation time to follow the temperature depen- 
dence predicted by the Landau -Teller theory at moderate 
temperatures. For temperatures in excess of about 8000K 
the vibrational model is extended to include an empirical re- 
sult for the relaxation time. The effect of introducing these 
temperature dependent collision numbers into the DSMC 
technique is assessed by making calculations representative 
of the stagnation streamline of a hypersonic space vehicle. 
Both thermal and chemical nonequilibrium effects are in- 
cluded while the flow conditions have been chosen such that 
ionization and radiation may be neglected. The introduc- 
tion of these new models is found to significandy affect the 
degree of thermal nonequilibrium observed in the fiowfield. 
Larger, and more widely ranging, differences in the results 
obtained with the different energy exchange probabilities 
are found when a significant amount of internal energy is 
included in the calculation of chemical nonequilibrium. 

Nomenclature 

AOTV = Aero-assisted Orbital Transfer Vehicle 
c = average molecular velocity 

C u = heat-transfer coefficient 

f(g) = relative velocity distribution function 
g = relative velocity of collision 

g* = characteristic velocity 

m r = reduced mass of collision 

n = number density 

T = translational temperature 

T* = characteristic temperature for Z R 

T re f = reference temperature for VHS model 

VHS = Variable Hard Sphere 
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X,* = mass fraction for species i 

Z R = rotational collision number 

(Z R ) oo = infinite temperature collision number 
Z t = translational collision number 

4>r = probability of rotational energy exchange 

<f> v = probability of vibrational energy exchange 
v = molecular collision rate 

u = parameter in VHS model 

cr re f = reference collision cross section for 
VHS model 

<j v = effective excitation cross section 
tit - Landau-Teller vibrational relaxation time 
Tp = Park’s empirical correction for relaxation time 

r v = total vibrational relaxation time 

Introduction 

The Direct Simulation Monte Carlo method (DSMC) de- 
veloped by Bird 1 is an important numerical technique for 
the calculation of low density flowfields. The method sim- 
ulates the flow at the molecular level by calculating col- 
lisions on a probabilistic basis. One of the most impor- 
tant aspects of DSMC calculations is the ability to model 
thermal nonequilibrium in which the temperatures asso- 
ciated with the various energy modes of the gas (trans- 
lational rotational, vibrational, electronic) are unequal. In 
order to reduce computational overheads, such phenom- 
ena are usually modelled in a simplistic manner in which 
the probability of energy exchange between these various 
modes is a constant for each type of transfer mechanism. If 
a particular collision is accepted for internal energy transfer, 
then post collision values are sampled from the local equi- 
librium distribution function in accordance with the Larsen- 
Borgnakke phenomenological model 2 . 

By employing this energy exchange scheme, significant 
degrees of thermal nonequilibrium have been calculated 
for expansions in plumes 3 ,4 and for compressions in shock 
waves 5 * 6 . Recently, one of the most interesting applications 
of the DSMC technique has been to the flowfield surround- 
ing a hypersonic space vehicle such as the Space Shuttle or 
an Aero-Assisted Orbital Transfer Vehicle (AOTV) as de- 
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scribed in Refs. 5 and 6. The degree of thermal nonequilib- 
rium in such energetic flows can have a significant impact 
on the shock standoff location, on the amount of chemi- 
cal activity, and on the convective and radiative heat load 
to the vehicle. It is therefore of great importance that the 
models employed in the calculation of energy transfer be 
physically realistic. For AOTV calculations the probabili- 
ties assumed for rotational and vibrational energy exchange 
are usually taken to be 0.2 and 0.02 respectively. The value 
for the transfer of rotational energy has been derived from 
experimental studies for temperatures of about 1000K. In 
the case of vibrational energy exchange, the probability of 
0.02 is an estimate of the maximum possible. In reality, 
the probability of energy transfer would be expected to be 
much smaller except at very high temperatures. It is the 
purpose of the present work to describe the procedures for 
implementing energy transfer probabilities which are de- 
pendent upon the local translational temperature. Rather 
than calculating the translational temperature at each cell in 
the simulation and then employing continuum expressions 
to evaluate the relaxation time 7 , each probability of energy 
exchange is obtained by integrating over all simulated col- 
lisions. This procedure should ensure that any nonequilib- 
rium attributes associated with the translational modes are 
incorporated into the calculation. The effect of introduc- 
ing such expressions into the DSMC technique is assessed 
by making one-dimensional calculations along the stagna- 
tion streamline of a hypersonic vehicle. First, the methods 
for calculating the probability of internal energy transfer are 
described. 

Probability of Rotational Energy Transfer 

The following approximate expression for the rotational 
collision number was obtained by Parker 8 
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where T* is the characteristic temperature of the intermolec- 
ular potential, and (Z^)^ is the limiting value. While 
Parker’s expression is derived from an analysis involving 
a large number of assumptions, the temperature dependent 
nature of Eq.(l) is in agreement with the more rigorous 
treatment of Lordi and Mates 9 who performed classical tra- 
jectory calculations. In the current work the value of (Zr)oo 
is chosen so as to obtain the best correspondence between 
Parker’s results and those of Lordi and Mates. To incorpo- 
rate Parker’s continuum formula into the DSMC technique 
an expression for the probability of energy transfer must 
be developed as a function of the relative collision velocity. 


This expression must reduce to Eq.(l) when integrated over 
all possible collisions. Such a formulation has been derived 
in Ref. 10 and is given as 
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where 4>r is the probability of energy transfer between the 
translational and rotational modes, Z t is the translational 
collision number, m T is the reduced mass of collision, and g 
is the relative collision velocity. In the derivation of Eq.(2) 
the Variable Hard Sphere collision model of Bird 11 has been 
employed in which the molecular collision rate is given by 

i/= 2no^ r(2 _ w) T]U (3) 
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where n is the number density, cr re / defines a collision cross 
section at a temperature T rc / y and u defines the interaction 
potential. In a separate set of calculations 10 , Eq.(2) was 
employed in the calculation of a standing shock wave and 
the results obtained were found to offer better correspon- 
dence to the experimental data than those with 4> R = 0.2. It 
has been found by Belikov et al 12 that Parker’s formulation 
appears to give the correct behavior up to about 4000K. In 
addition, it has been shown by Lumpkin et al 13 that Eq.(l) 
gives good agreement with experimental shock wave thick- 
ness for a range of Mach number. Indeed, in Ref. 13 it 
is observed that Parker’s simple formulation provides bet- 
ter correspondence than the results obtained with calcula- 
tions which employed a more sophisticated model based 
on rotational transition rates. As shown in Fig. 1, even at 
1000K the probability of rotational energy transfer for ni- 
trogen is certainly different to 0.2 and the variation in <j) R 
will therefore have some impact on the calculated results. 
When Eq.(2) is evaluated at each collision, it is found that 
equipartition of the thermal modes is not achieved due to 
the fact that the collision probability is biased towards those 
collisions having smaller translational energies. The mean 
collision probability for each cell in the computational do- 
main is therefore employed, and is obtained by averaging 
Eq.(2) over all collisions. 

Probability of Vibrational Energy Transfer 
In a pure gas of diatomic molecules, the mean proba- 
bility of vibrational-translational energy exchange may be 
expressed as 

<t>v = — = [ 4> v (g)f(g)dg (4) 

t v v 



where 4> v is the probability of energy exchange for the col- 
lision with relative velocity g, f(g) is the relative velocity 
distribution function, and r v is the vibrational relaxation 
time. In the following, we wish to choose a form for $ v 
such that the vibrational relaxation time is given by the ex- 
pression obtained by Millikan and White 14 who correlated 
a number of experimental results. The probability of energy 
exchange is assumed to be of the form 

<f>v = -l-g a exp( — — ) (5) 
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where a and Z 0 are constants, and g* is a characteristic ve- 
locity. The form of Eq.(5) is identical to that originally as- 
sumed by Landau-Teller with a set equal to zero. The in- 
tegral obtained by substituting Eq.(5) into Eq.(4) can only 
be evaluated approximately using the method of steepest 
descent 15 . The result obtained with this method reveals that 
the temperature dependent nature of the vibrational colli- 
sion number is obtained when 

a = 3 — 2uj (6) 


A detailed description of the preceding analysis together 
with the expressions for the constant Z 0 in Eq.(5) and the 
characteristic velocity g* are included in Ref. 15. For a 
hard sphere interaction (w=0), the instantaneous probability 
of vibrational-translational energy exchange for nitrogen is 

4^00 

<f> v = 1.20 x 10 ~ n g 3 exp( ) (7) 
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Results obtained with the DSMC technique employing 
Eq.(7) are plotted in Fig. 2 together with the continuum 
result. Due to the approximation used in evaluating Eq.(4), 
the results are not identical, with the DSMC results always 
lying less than a factor of 2 above the exact solution. 

The vibrational relaxation time of Millikan and White is 
generally recognized to be valid up to temperatures of about 
8000K. In the consideration of hypersonic flow surrounding 
space vehicles, the translational temperature may be as high 
as 40,000K. To better simulate the vibrational-translational 
exchange process at such elevated temperatures Park 16 in- 
troduced the following empirical correction 

T v = T LT + T p (8) 

where tl T is the Landau-Teller relaxation time. 


a v is the effective excitation cross-section and c is the av- 
erage molecular velocity. This correction to the vibrational 
collision number may be evaluated in a similar manner to 
that employed for rotational-translational energy exchange. 
If the excitation cross-section is taken to be constant, then 
the following expression is obtained 

^ v£.,r(2- u ) [r / (2 -“ )r "' 1 “ < 10 > 

An expression which reduces to Eq.(10) when integrated 
over all collisions is 

6 = ° v \ mrg2 
P y/la nJ 1.2(2 -uj)kT Tef 

It should be noted that it is possible to include a tempera- 
ture dependent form for cr v if desired. From Eq.(ll) it is 
seen that for hard spheres, the result is independent of the 
collision velocity. Park’s modification to the vibrational re- 
laxation time is implemented in the DSMC method by av- 
eraging Eq.(ll) over all collisions. To ensure that the inte- 
gration is being performed correctly, the results for u»=0.25 
and nitrogen molecules are shown in Fig. 3 for the modi- 
fied vibrational relaxation time of Eq.(8). It is seen that the 
DSMC and continuum results are in excellent agreement. 

Calculations 

Having derived expressions which allow both the rota- 
tional and vibrational collision numbers to vary with tem- 
perature, the effect of introducing these developments into 
energetic DSMC calculations is now assessed. The prob- 
lem chosen for consideration is the flow along the stagna- 
tion streamline in front of a hypersonic space vehicle. The 
solution procedure follows that of Bird 17 who has shown 
that the stagnation streamline flow may be reduced to one 
dimension. One end of the computational domain is open 
through which enter molecules representative of the up- 
stream conditions. At the opposite end is the surface of 
the space vehicle, which is assumed to be a diffuse reflec- 
tor with full energy accomodation. As the molecules are 
reflected back from the wall a shock wave forms which 
gradually moves away from the surface. A steady flow is 
achieved by removing molecules from a downstream por- 
tion of the flow with probability proportional to the square 
of the normal velocity component. This procedure has 
been shown to conserve energy and momentum when the 
molecules are removed in the region lying between the 
shock standoff location and the surface. 

The conditions investigated represent the exit of the 
space vehicle through the atmosphere rather than the reen- 
try. The exit phase is significantly less energetic than reen- 
try and allows us to neglect ionization and radiation effects 




in the calculations. An altitude of 90km has been selected 
and the free stream flow conditions are listed in Table 1. 
The wall surface is assumed to be constant at 1000K and is 
non-catalytic. A five species gas mixture is employed (N 2 , 
O 2 , NO, N, O) with 19 different chemical reactions (see 
Table 2(a), Ref. 18). Due to the low densities involved, 
recombination reactions have been omitted from the anal- 
ysis. Calculations have been completed for three different 
cases. The first of these follows the previous work by set- 
ting 0.2 and <£y=0.02. Secondly, the temperature and 
species dependent expression in Eq.(2) for rotational en- 
ergy transfer is introduced. Finally, both the rotational and 
the vibrational energy exchange probabilities are calculated 
as a function of temperature and molecular species using 
Eq.(2) and Eqs.(6,ll). Details of the parameters used in the 
energy exchange models are given in Table 2 for rotation, 
and in Table 3 for vibration. The value for the translational 
relaxation number Z t> required in Eq.(2), is that proposed 
by Fritzsche and Cukrowski 19 . 

While the primary focus of the present work is the inves- 
tigation of thermal nonequilibrium effects, the importance 
of chemical nonequilibrium cannot be ignored. Initially, ro- 
tational energy contributions to the energy available for re- 
action are only included on a restricted basis, and where re- 
quired by Bird’s steric factor for the VHS model. The frac- 
tions of rotational energy allowed to participate in the var- 
ious reactions are those employed by Bird 20 in his general 
purpose code. The effect of this assumption on the flow- 
field is also analysed by allowing all of the rotational energy 
available in each collision to be incorporated into the steric 
factor. For the dissociation reactions it would of course be 
more realistic to include contributions from the vibrational 
modes. Unfortunately, the implementation of such a pro- 
cedure is numerically expensive when Bird’s steric factor 
is employed. It is proposed that the use of the rotational 
energy will at least give an indication of the importance of 
including internal energy contributions in the calculation of 
chemical nonequilibrium. 

Results and Discussion 

At an altitude of 90km, Ref. 16 shows that the shock 
standoff distance, which is defined as the point at which the 
local density is 6 times the free stream value, is located at 
a distance of 0.11m from the vehicle surface. The present 
calculations have been performed by forcing the standoff 
distance to be 0. 1 1m for the case where the exchange prob- 
abilities are constant The other cases are then run under 
identical conditions such that the steady state is introduced 
at the same point in the simulation. This allows investiga- 


tion of any movement in the standoff distance in the various 
calculations undertaken. In Fig. 4a, the temperature pro- 
files along the stagnation streamline are plotted for the case 
in which the probabilities of rotational and vibrational en- 
ergy exchange are both kept constant. As expected, a signif- 
icant degree of nonequilibrium exists, with the translational 
mode reaching a peak temperature of about 25,000K. The 
results shown offer excellent agreement with the profiles re- 
ported in Refs. 5 and 6. The variation along the stagnation 
streamline of the mass fraction of each of the five species is 
shown in Fig. 4b for the same case. A significant amount 
of chemical activity is observed in the flowfield. 

Let us now consider the effect of introducing a tem- 
perature and species dependent rotational energy exchange 
probability into the calculations. In Fig. 5 the rotational 
temperature profiles obtained with the two different mod- 
elling methods are shown. It is found that the maximum 
rotational temperature attained is significantly reduced for 
the case of the variable exchange probability. This effect 
may be explained by examination of Fig. 6 in which the ro- 
tational probabilities obtained in the calculations are shown 
for each of the three molecular species. It may be seen from 
Fig. 4b that molecular nitrogen dominates the flow, so let us 
concentrate on this data. At the relatively low free-stream 
temperature, the probability of rotational energy transfer is 
very close to the nominal value of 0.2 and gradually de- 
creases as the translational temperature increases. It should 
be noted that the energy exchange probability increases 
again as the cooling effects of the vehicle surface become 
important The values of the exchange probabilities for 
oxygen are observed to be everywhere greater than those for 
nitrogen with a maximum of about 0.29 at the free-stream 
boundary. For nitric oxide, a value for the exchange prob- 
ability is not obtained until a distance of about 0.2 from 
the vehicle. This is explained by the fact that, up to this 
point, NO gas exists in insufficient quantities to produce any 
collisions in the simulation. As the rotational energy ex- 
change probability is quite large at the upstream boundary 
for both the nitrogen and oxygen molecules then it might be 
expected that the rotational temperature in Fig. 4a would lie 
closer to the translational result at this location. However, 
Fig. 7 reveals that the various thermal modes are in fact in 
equilibrium at the upstream location for nitrogen, and this 
is also found to be the case for the oxygen molecules. Ex- 
amination of the results obtained for each species revealed 
that the large translational temperature obtained in Fig. 4a 
at the upstream boundary is caused by a very small fraction 
of highly energetic, backscattered atoms and nitric oxide 
molecules. 



In Fig. 8, the vibrational temperature profiles are shown 
for the two energy exchange probability models. As in the 
case of rotational relaxation, it is found that the maximum 
vibrational temperature is significantly reduced. The in- 
dividual exchange probabilities for each molecular species 
are plotted in Fig. 9. Careful attention should be given to 
the labelling of the axes. Only a portion of the domain mod- 
elled along the stangation streamline is included. Beyond 
distances of about 0.25m from the vehicle, the collision rate 
is too low to allow any vibrational energy exchange for 
the variable probability model. It is to be remembered that 
the DSMC technique employs a discrete number of simu- 
lated molecules in the simulation which undergo a discrete 
number of collisions. For calculated exhange probabili- 
ties which are very small, e.g. 10“ 4 , it is therefore diffi- 
cult to accurately resolve the exchange phenomenon. In di- 
rect contrast to the results found for the rotational exchange 
probabilities, the maximum translational temperature gives 
rise to the maximum exchange probabilities. For both ni- 
trogen and oxygen these maximum values fall well below 
0.02 which is the value previously assumed. It is satisfying 
to note that the values for nitric oxide lie well above those 
of the other molecules. This is the trend expected from ex- 
periment, and shows that such behavior can now be incor- 
porated into the DSMC calculations through application of 
the variable exchange probability model. 

It is therefore found that the introduction of the variable 
exchange probabilities leads to a greater degree of thermal 
nonequilibrium in the flow due to the fact that the calcu- 
lated probabilities are always much smaller than the values 
of 0.2 and 0.02 normally assumed for energy exchange with 
the rotational and vibrational modes respectively. However, 
at the macroscopic level, only a small effect is observed 
through the introduction of the temperature dependent mod- 
els. One example is shown in Fig. 10 where the mass frac- 
tion of molecular and atomic nitrogen is compared. It is 
seen that a greater amount of nitrogen atoms are produced 
in the case where variable probabilities are employed. In 
Table 4 various results obtained from these simulations are 
shown. It is found that both the shock standoff distance and 
the heat coefficient increase slightly with the introduction of 
first the variable rotational exchange probability, and then 
additionally for the corresponding model for vibration. 

It is perhaps surprising that the introduction of the new 
models has had such little effect on the important flow quan- 
tities. However, this is almost certainly due to the fact that 
the internal energy modes are only included in the chemi- 
cal nonequilibrium phenomenon on a restricted basis. At a 
time when the significant role of internal energy in chem- 


ical nonequilbirium has been incorporated into continuum 
calculations 21 <22 it is unsatisfactory for the discrete particle 
solution methods to lag behind in this area. There is a real 
requirement for improvement in the modelling of chemical 
reactions in the DSMC technique. In an attempt to assess 
the effect of the particular chemical model employed, the 
same flow conditions have been recalculated with the total 
rotational energy now being available for reaction. 

In Fig. 11, the mass fraction of molecular nitrogen is 
shown for the cases where constant exchange probabilities 
were employed for both rotational and vibrational energy. 
It is seen that a reduced amount of dissociation is found for 
the increased internal energy contribution. This would be 
expected as Bird’s steric factors assume perfect equilibrium 
for the various thermal modes. Therefore, when the internal 
temperatures lie below the translational result, which is the 
case for the flow along most of the stagnation streamline, a 
smaller rate of reaction will occur. 

The thermal nonequilibrium effects discussed above for 
the restricted rotational energy contribution to chemical re- 
actions were repeated in the new set of calculations. How- 
ever, larger differences for other flow quantities were dis- 
cernible. In Fig. 12 are shown the mass fractions for atomic 
and molecular oxygen for the two cases where constant and 
variable exchange probabilities were employed. There is a 
smaller amount of molecular oxygen present in the mixture 
when the variable models are introduced. Similar behavior 
is observed for nitrogen, as shown in Fig. 13, but to a lesser 
extent The exact explanation for these trends is compli- 
cated by the presence of the number of competing reactions 
which are accounted for in the calculations. The differences 
between the solutions obtained with the variable exchange 
probabilities are now large enough to be noticed in such 
macroscopic quantities as density and velocity which are 
plotted in Figs. 14 and 15 respectively. It is found that the 
density is greater for the variable probability calculations 
in the region lying between the free stream and the shock 
standoff location. Thereafter, the variable result is slightly 
below that obtained with the constant probability. In Fig. 
15, it is seen that the velocity solutions obtained with the 
different methods are clearly different with the variable re- 
sult lagging behind. These findings are of a similar mag- 
nitude to the differences noted by Moss et al in Ref. 18 
for different sets of chemical reaction rate constants. It is 
therefore found that further uncertainty in such calculations 
arises from the techniques employed in the simulation of 
thermal and chemical nonequilibrium effects. 

Various global results obtained with these calculations 



are compared in Table 5 with that previously obtained for 
the constant exchange probabilities. The assumption of a 
full rotational energy contribution for chemical reactions re- 
sults in an increased effect through the inclusion of the vari- 
able energy exchange probability models. The heat-transfer 
coefficient, the mass fraction of molecular nitrogen at the 
surface, and the shock standoff distance are all affected for 
the different modelling methods. 

Concluding Remarks 

The calculations reported show that the introduction 
of temperature and species dependent internal energy ex- 
change probabilities into DSMC calculations can have sig- 
nificant impact on the results obtained. The constant val- 
ues normally employed for rotational and vibrational en- 
ergy transfer are generally larger than the values calculated 
in the simulation. The degree of thermal nonequilibrium 
is therefore increased due to a reduction in the amount en- 
ergy transfer performed. Indeed, the vibrational exchange 
probabilities obtained for nitrogen are an order of magni- 
tude smaller than the constant value previously assumed. 

The differences observed in other flow quantities for the 
different energy exchange models are found to be largest 
when the internal modes are allowed to contribute a sig- 
nificant amount of energy to that available for each chemi- 
cal reaction. In the current investigation only the rotational 
energy contribution has been considered, and is found to 
be important. In the case where the total rotational energy 
is included in the collision energy, the variable energy ex- 
change models are found to have a pronounced effect on 
such important quantities as the shock standoff distance and 
the heat-transfer coefficient. 

It is concluded that consideration should be given to 
the implementation of the energy transfer probabilities de- 
scribed in this paper. It is probable that the effects observed 
in the present calculations will be more pronounced for the 
more energetic flowfields associated with atmospheric re- 
entry of an AOTV or Mars return vehicles. It is a require- 
ment for the future to investigate such conditions and to 
consider in detail the role of chemical nonequilibrium in 
such flows. In addition, the role of the variable exchange 
probabilities in expanding flows, such as that around the 
aerobrake skirt of an AOTV, should be considered. In such 
flow, the freezing of the vibrational temperature plays a sig- 
nificant role in the determination of radiation effects. 
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Altitude (km) 


Uco (km/s) 


Table 1. Freestream conditions. 

Density (kg/m 3 ) (K) 




Xo 2 


90 

7.5 3.418xl0- 6 

188.0 

0.23 

0.77 

Table 2. Parameters for Rotational Energy Exchange. 

Species 

(ZR) oo 

T* (K) 

-Zt 


n 2 

23.3 

91.5 

1.5 


0 2 

16.5 

113.5 

1.5 


NO 

19.5 

100.0 

1.5 



Table 3. Parameters for Vibrational Energy Exchange. 



Species 

z 0 

g* (m/s) 

a„ (m 2 ) 


N 2 

8.30x10'° 

43300 

10-2° 


o 2 

1.50xl0 n 

18200 

S 

1 

o 


NO 

1.51xl0 u 

6400 

S 

i 

o 

Table 4. Results For Restricted Rotational Contributions To Chemical Reactions. 

4 >r 

4>V 

Shock Standoff (m) 

C H 

Surface Mass Fraction of N 2 

0.2 

0.02 

0.110 

0.177 

0.561 

Eq.(2) 

0.02 

0.113 

0.178 

0.547 

Eq.(2) 

Eq.(8) 

0.115 

0.180 

0.540 


Table 5. Results For Different Rotational Contributions To Chemical Reactions. 


Rotational Contribution 

4 >r 

4>v 

Shock Standoff (m) 

C H 

Surface Mass Fraction of N 2 

Restricted 

0.2 

0.02 

0.110 

0.177 

0.561 

Full 

0.2 

0.02 

0.110 

0.177 

0.552 

Full 

Eq.(2) 

0.02 

0.116 

0.182 

0.530 

Full 

Eq.(2) 

Eq.(8) 

0.118 

0.188 

0.506 







Temperature (K) 


Fig. 1. Rotational energy exchange probability as a func- 
tion of temperature. 



Fig. 2. Vibrational energy exchange probability as a func- 
tion of temperature. 



Temperature (K) xIO 4 

Fig. 3. Vibrational energy exchange probability with Parks 
empirical correction. 


xIO 4 3.0 
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Fig. 4a. Thermal nonequilibrium along stagnation stream- 
line for constant exchange probabilities. 



Fig. 4b. Species mass fractions along stagnation streamline 
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Fig. 5. Comparison of rotational temperature profiles for 
different rotational energy exchange models. 
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Fig. 7. Thermal nonequilibrium of molecular nitrogen. 



Fig. 8. Comparison of vibrational temperature profiles for 
different vibrational energy exchange models. 
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Fig. 10. Comparison of mass fractions of molecular and 
atomic nitrogen for different energy exchange models. 



Fig. 11. Comparison of mass fraction of molecular nitrogen 
for different stcric factors. 
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Fig. 12. Comparison of mass fractions of molecular and 
atomic oxygen for different energy exchange models with 
full rotational contribution to the steric factor. 


Fig. 14. Comparison of density for different energy ex- 
change models with full rotational contribution to the steric 
factor. 



Distance (m) 

Fig. 13. Comparison of mass fractions of molecular and 
atomic nitrogen for different energy exchange models with 
full rotational contribution to the steric factor. 
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Fig. 15. Comparison of velocity for different energy ex- 
change models with full rotational contribution to the steric 
factor. 




